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Estimators of outer scales of atmospheric turbulence usually fit the phase screen snapshots derived 
from local wave front sensors to a Zernike basis, and then compare the spectrum of expansion 
coefficients in this basis with a narrowing associated with decreasing outer scales. 

This manuscript discusses aspects that arise if the Zernike basis is exchanged for a Karhunen- 
Loeve basis of statistically independent modes. Data acquisition turns out to be more demanding 
because sensing the tip-tilt mode is required. The data reduction methodology is replaced by fitting 
of variance ratios of an entire (long exposure) data set. Statistical testing of hypotheses on outer 
scale models can be applied to a set of modes — supposing other noise originating from detector 
readout and the optical train can be disentangled. 
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I. STATISTICS OF ATMOSPHERIC TURBULENCE 

^ r"| A. Measuring and Interpreting Pupil Phases 

A contemporary adaptive optics (AO) system of an astronomic telescope measures angles of arrival of the incoming 
wave front sampled across segmented areas that subdivide the full entrance pupil. The task is to weave the two 
principal parameters — motion of the image either represented by two Cartesian directions on the detector or by an 
angle of arrival and azimuth [23| — of the sensor array into a phase function covering the entire pupil. For the purpose 
of the present work we leave aside aspects of phase unwrapping [l9l l29l [33j , finite element discretization 34] , or noise 
CsJ ■ besides the atmospheric turbulence itself 

As the information obtained represents the phase in the pupil, optical path differences integrated along paths 
through the entire atmosphere and altitudes dependencies are basically not available within a single AO system. 
Thicknesses and altitude variation of the structure constants may vary Understanding of the origin of finite 

outer scales of turbulence remains poor for that reason, and the scattering of outer scale lengths at wavelengths in 
the visible is significant [3, Ilia, [la |38| . 

Optical path differences from stellar interferometers with baseline lengths of the order of the outer scales promise 
higher sensitivity. If the interferometric subsystems do not share any data streams with the AO systems of the 
individual telescopes, the time series of the fringe tracker may be the only information available, and then mediation 
by some wind velocity with its own set of unknowns is needed to fit the shape of temporal spectra @, [13] to outer 
scales. Results from the Palomar Testbed Interferometer merged with tip-tilt information of the telescopes show outer 
scales of some tens of meters — compatible with data from single-telescope apertures [28| . 

Auxiliary data products of this kind are not a standard of reduction pipelines of interferometric observations for 
now, so we discuss only the case of a single circular aperture. The remainder of Section Q] summarizes the literature 
on the von Karman model of outer scales with a single layer with isotropic turbulence. Section [II] introduces its 
modal composition for the circular pupil. On the technical side, the matrix elements that arise in the Zernike basis 
are reduced to Hypergeometric Functions (Section III A|) . and central obscurations common to on-axis mirror trains 
are handled in a hybrid form that separates Fourier and real space with the aid of a Neuman series (Section III C[) . 
Benefits and problems compared to the use of Zernike covariances are shortly discussed in Section Hi Bl 
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B. The von Karman Model of Outer Scales 

The refractive index structure function T> n as a function of the 3D separation Ar, and its Fourier representation 
$„ at spatial frequencies / are related via 



O n (Ar) = 2 J $„(/)[! - cos(27rf • Ar)]d 3 /- 



The von Karman model of the power spectrum in particular is 

*«(/) = c„C 2 (/ 2 + / T (7+3)/2 



(1) 



(2) 



where C 2 is the structure constant which characterizes the global strength of the turbulence, and where fo is a 
wavenumber equal to the inverse of the outer scale. 7 is a power spectral index which for our analysis will be set to 
the Kolmogorov value of 2/3; formulas will generally not be reduced with this value to support independent tests of 
this scaling [|, 0, El, H3, E3| • The constant 



27r 3/2+ 7r (^7) 



0.00969315 (7 = 2/3) 



is settled by demanding that at small Ar 



P„(Ar) = C 2 n \Ar\ 



(3) 



(4) 



This is a model of isotropic turbulence, so the information can be bound to the radial dependencies in real and 
Fourier space by integration over polar and azimuth angles, 



P„(r) = 8^ ^ / 2 [1 - j (2nfr)]$ n (f)df = 87rc„C 2 ^ / 2 [1 - j (2wfr)} (/2 - /2)(3+7)/2 d/- 

Performing the integration leads to a representation in terms of Modified Bessel Functions K [111, [l2|, [H| , 

2 



± V Wr -7 - (Wo)^r(7/2) 

cr n(r)r r(- 7 /2) 



1 



r(7/2) 

The limit of infinite outer scale recovers the Kolmogorov power law 

which is the value reached by the curve of r = R in the left graph of Fig. [T] 



(irrf Q y /2 K l/2 (2irrf ) 



(5) 



(6) 



(7) 



Phase Structure Functions 



In most applications, the observable is not the three-dimensional distribution of the refractive indices but the phase 
difference tp between light rays at an optical wavenumber k (which is 2tt divided by the wavelength) . We write T> v for 
the structure function of the phase if observed over a projected distance P in an input pupil of the instrument after 
two rays have traveled along parallel lines through a layer of height h (which includes the air mass factor) in which 
the refractive index structure function modeled in the previous section defines the statistics. This type of projection 
onto a two-dimensional plane introduces a weighting of the refractive index power spectrum <£>„ with Bessel Functions 



V V {P) = AnFh / [1 - J {2nfP)]$ n (f)fdf, 



for the von Karman spectrum ^ condensed into the well-known 



V V {P) = -k 2 hClP l+ ~i 



r(-7/2)( 7 rP/ ) 1 



1 + 7 



- 2(7rP/ )( 1+ ^/ 2 ^ (1+7)/2 (2^P/ ) 



(8) 



(9) 
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FIG. 1: The refractive index structure function T> n (r) / 'C„ of ([6]) at 7 = 2/3 as a function of outer scale with the radial distance 
r as the parameter, or with the roles of these two dependencies swapped. T> n depends essentially on the product rfo, so scaling 
with some anonymous reference length R leads to the universal curves shown. 



The Kolmogorov limit of infinite outer scale is the familiar 



where 



7(7) 



lim V V {P) = PhCftP^gfr) 

JO— HJ 



'7f r( J/ 2 ,1 ,2) « 2.9143808, (7 = 2/3). 



iX-7/2) 



(10) 



(11) 



Collecting variances of optical path differences from this model — obtained by striking the factor k 2 of V v — we 
notice that the strength of the fluctuations remains smaller than predicted by the factor 3(7) in (fTTj) . because this 
has been derived in the limit h/P — > 00 of infinite layer height [24|. For finite h/P, 



2k' 



(h-Az) V n {^{Az) 2 + P 2 )-V n {Az) d(Az) 



2k 2 ClP' t+1 h 



h/P 



h/P 



(t 2 + \y' 2 - f 



dt. 



where the integral can be rephrased as [2^, (3. 243.1), (3. 251.1)] 



5(7) = 2 
with [H, Hi 



h/P 



[l-t/(h/P)] (1 + t 2 )^/ 2 -V 



dt 



P 



3F2 



-7/2,1/2,1 
3/2,2 



3F2 



-7/2,1/2,1 
3/2,2 



z) = 2 2 F, 



-7/2,1/2 
3/2 



-7/2,1 

2 



(h/P) 



1+7 



( 7 + l)(l+7/2)' 



(12) 
(13) 

(14) 

(15) 



This factor cannot reach 2.91 in practise: the layer height h multiplied by the air mass is limited to the order of 10 km 
representing the entire atmosphere, and P is of the order of 2 m for a single telescope up 100 m for an interferometer. 
Then h/P is in the range 5000 to 100, and Figure [5] demonstrates that this finite layer thickness places 5(7) in the 
range 2.2 to 2.6. Long baseline interferometers "miss" some power of T> v through this geometric sampling effect. 
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FIG. 2: The undersampling factor (|14l) as a function of the ratio h/P. 



II. KARHUNEN-LOEVE MODES IN THE ENTRANCE PUPIL 

A. Phase Structure Function in the Circular Pupil 

Removal of the spatial mean of the structure function ([5]) defines associated covariances C v for two-dimensional 
vectors r in the pupil plane and two-dimensional spatial wavenumbers /, 



where 



C v (t) = I C v (f)e~ 2wlfr d 2 f = 2tt J C v (f)J (2wfr)fdf, 



C 9 (f) = k 2 h$ n (f). 



(16) 



(17) 



The statistically independent Karhunen-Loeve (KL) modes Fj are eigenvectors of the integral operator [2(| 

/ C v {v-v')F j {v')d 2 r' = B' i F j {v), (18) 

J\r'\<R 

where the integration is over the area limited by the pupil radius R. Moving on to Fourier space, the convolution 
simplifies to 



C v (J)F j {f)=B i F j {r) i 



and more explicitly with (1171) to 



n "(/ 2 + /o 2 ) ( ^ +3)/2 



F J (f)=S 2 ^(f) 



(19) 



(20) 



Introducing dimensionless variables r = i?x and i?f = er helps to gather a set of common constants to unify the 
analysis, 



1 



(ff 2 +0-2)(T+3)/2 



B 2 



k 2 hC?,Rf+ 3 



The notation assumes that azimuthal dependencies have been split off, 



(21) 



(22) 



Mq {9) = M x ( C0S ^ ; e , s ( J' f = °' 
9W V 2tt 1 sm(q9) q \ 1, \q\ > 1, 



(23) 
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such that the two discrete indices p and q can be used to enumerate the modes by number of nodes in the radial 
direction and periodicity along the azimuth 9. If the Kp(x) are expanded in the basis of Zernike polynomials i?^, 

K i q \x)= E r n , p , q y/2(n+l)Rl(x), 

n=q mod 2 

Noll's Fourier representation of the induces a Fourier representation of the radial component [32l l36j, 
K^(a)= £ r», p , g V2^TT)(-l) ( "^ /2 » ^" +l(2 " a) = £ ^.V^+l)*" 

' fT ' (7 

n=g mod 2 n=g mod 2 

which becomes an eigenvalue equation for the vector of the expansion coefficients T n ,p,q (3lj |. 

dk B 2 



(24) 
(25) 

(26) 



P CO 

c„(27r)T+ 3 V r n , p , gV /2(^TT)i"- n ' / J n+1 (k)J n , +1 (k) 



k(k 2 + k 2 )^)/ 2 k 2 hC^R-r+ 3 '' n '' p ' 9 'y/2{n' + 1) 



n— g mod 2 



At this stage some arbitrariness remains in scaling both sides of the equation, depending on preferences of normal- 
ization to pupil areas, squared diameters or radii. I shall divide through 3(7) to gather 5(7) in front of the C 2 of the 
right hand side to match (fTOj) . and shall multiply by a factor 



2c^ 2 



4r 



3 + 7 
1 + 7 



(l+7)/2 



24 

y r(6/5) 



5/6 



6.883877, (7 = 2/3) 



for compatibility with calculations that introduce a Fried radius as a reference. Altogether, 

r(i + f) 



2 2 t n=a mod 2 



defining 



and 



2 _ 2c y B 2 
P ' q ~ 5(7) k 2 hCl{2R)i+* 



Jn+l{k)- 



k(k 2 + fe2)(3+7)/2 

This integral is evaluated as proposed by Winker [38|, , 



Inn' — 



dk 







tt/2 



d9J n + n '+2 (2fc cos 0) cos[(n — 



J n , +1 (k)dk. 



1 







00 /' 7r /2 1 

dfc / rf6> — 
27ri 



2tTZ 



fc(fc 2 + fc 2 )( 3 +^)/ 2 

r(-s)(fccos6>)"+"' +2+2s 

— — L — : cos ti — n )9\ 

7iT(n + n' + 4 + 2s) 



k(k 2 + fc 2 )( 3 +^)/ 2 



(28) 
(29) 

(30) 
(31) 

(32) 
(33) 



ioo Y( — S ^k n+ +2+2s 

_ loo dS T(n + ri + 3 + s) 2 n + n '+ 3 + 2s {n + ri + 3 + 2s)T{n + 2 + s)T(n' + 2 + s) k(k 2 + k%)^+i)/z 



„n+n — 1— 7 



r ( l+7-"-n' ) r (n±nl + ^ 

2™+"'+ 3 r(n + 2)r(n' + 2)r(^p) 



+ 



„ . n + ri n + ri n + ri + 3 , , 

x 3^4(— ^— + 1, — ^— + 2, ; n + 2, n' + 2, n + n' + 3, 

r( "+"' 2 - 1 ~ 7 )r(4 + 7) 



l + n + n'- 7 2 



Q"y+4p^ 7+7+w+n' ^p^- 5+f+n-n' 5+7+w' — 71 ^ 



,5 + 7 3 + 7 7 5 + 7 + n — n' 5 + 7 + n' — n 7 + j + n + ri 3 + 7 - n - n' 2 
x 3 ^4(^-,^-,2 + -; - , - , - , - ;fc ) 



(34) 



6 



Prototypical results of the calculation are listed in App. [AJ The reduction of the eigenvalues as the outer scale 
shrinks is visualized in Fig. [3J and the inward motion of the normalized radial functions Kp (x) is shown in Fig. |4l 
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0.001 




FIG. 3: The 
< a < 0.8. 



K < 1 > K < 2 > K <°> K < 3 > K/ 1 » K < 4 > K/ 2 > k/°> K < 5 > k/ 3 > K 2 < 1 > K < 6 > 
eigenvalues Xp >g (|30[l of the KL modes of the phase variance as a function of the outer scale wavenumber 



The modes at the largest scales are effected most, as known from earlier analyses |4l| 



B. Estimation of Outer Scales 



The customary reduction of AO data to outer scales reconstructs the phase as a function of the location in the 
pupil plane, collects the expansion coefficients of the Zernike expansion, and compares the relative strengths of the 
coefficients with the reduction predicted as a function of the outer scale [U [41| . 

Which new aspects arise if a chart of the KL modes as in Figure [3] is used instead? 

Obtaining the expansion coefficients still is a matter of fitting a phase function to a set of modes. One may start with 
a Zernike basis [3(| and convert its coefficients to the KL basis with the aid of Appendix [A] The appendix provides 
the matrix that converts from the Zernike to the KL basis, but no additional computation is demanded because the 
inverse matrix is the transpose matrix. At this step, the tip-tilt component is an essential ingredient to the process — 
this is a requirement that was not known to the approach based on the Zernike covariances, and probably disables 
processing if tip-tilt stages are autonomous and data not available. If a modular approach to information processing 
of various layers of phase corrections in telescope operation — from course with fringe tracking over medium of tip-tilt 
stages to fine of adaptive optics — results in a lean system optimized to dispose intermediate data streams irrelevant 
to the observing astronomer, post-processing of the engineering data concerning components like the atmosphere is 
put on a diet. 

Another modest technical disadvantage is that fitting process becomes iterative in nature because the entries of 
Appendix |A1 depend on the outer scale, which is the unknown. This is only a matter of producing this static table for 
a range of outer scales of interest and therefore no real concern. 

The conceptual difference is that the data reduction is based on the statistical independence of the KL modes. The 
outer scale is selected that generates measured variances of the modal coefficients proportional to their Xp q . (The 
overall factor C\ drops out of the calculation.) The technique does not derive the outer scale from any individual 
snapshot, but it starts from a hypothetical outer scale, reduces the series of exposures to variances of modal coefficients, 
and readjusts the value of the outer scale until these variances match that common statistics. In addition, standard 
statistical methodologies can be applied to derive confidence intervals of the result. 
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FIG. 4: Influence of the inverse outer scale < ao < 1 on some large-scale KL modes (|25|) of the phase covariance in the 
circular input pupil. 



C. Central Obscuration 



A central obscuration from a secondary mirror removes data up to a radius e [measured in units of the pupil radius, 
< e < 1, unrelated to the e of (|23p ] from the areal statistics. The radial component of (p?Hj) is represented by the 
Fourier pair 

/l />oo 
xdxK( q \x)J q {2irax), K^\x) = 2i:{~i) q J uduKf (a)J q (27rax), (35) 
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which recovers the results of the previous section if e = 0. Insertion of Kp q \cr) in (|2"Tj) yields the eigenvalue equation 
in mixed coordinates, 

f°° 1 B 2 

M-i)^ I (ff2 + <7o2)(7+3)/2 ^ ) WJ 9 (2^)da = j^^K^x). (36) 

Insertion of (|35[) and scaling the wave number with k = 2ira means 

00 i r 1 , s B 2 

x 111 t-r.r(n\t f\i t ft t\ t ft \ *-* 



W +3 c n ] o ^ (fc2 + fco2)(7+ 3)/ 2 / ^(W,W = ^^). (37) 
The integral over k is managed by isolating variables 0[3 (§7-15)] 

\j q {kx)J q (kx') = {xx'Y £)(-l)»(2« + 2n + l)J 2 , +2n+ i(fc) ^' " ^ / "■ " 



n>0 
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/ -ra, 2g + n + l I . 2 i\ tp ( -n,2q + n + l , . 2 ,\ /ooN 

x 2*1 I 9 + j | sm \ 2 Fi I ^ + x | sin $ I . (38) 

Here, the two auxiliary variables </> and $ are defined by 

cos cos $ = x; sin ^ sin $ = x', (39) 

that is 



sin20 = l-{x 2 -x' 2 ) + ^l-{x + x') 2 ^l-{x-x') 2 ^ 



,^ 1-(X 2 - X 12 ) - y/1 - (X + X') 2 y/1 -(X- X'f 



sin 2 $ = - — — . (41) 

The two Hypergeometric Functions in (|38]l are Jacobi Polynomials [1, (22.5.42)]. Each term is then integrated with 
[H 6.565.8] 

J 2g+2 »+i(fc) ,, _ ^^(H^^f-g-n) / l + q + n fc| 

o (k 2 + k 2 )(^)/ 2 4 9+«+ir(2±3)r(2g + 2n + 2) 1 2 \q + n + 2n + 2q + 2 1 4 



To avoid that x + x' > 1 in (|4Tj) . another variable transform t = 2fc is added. The following results have been obtained 
with a finite element decomposition of the interval [e, 1], and introducing the function 

KW(x) = ^/2K^(x), (43) 
such that the matrix on the left hand side, which couples each x to all the x', becomes symmetric: 

The piston mode has been reduced as outlined in Appendix [Bj 

If the radius is selected equivalent to the M2/M1 ratio of the telescopes of the Very Large Telescope, Figure [S] 
results. Besides the absence of values in the range x < e, these shapes are hardly different from the full aperture 
modes of Figure |4j 

The variances A 2 in the individual modes are shown in Figure [5] [They are still defined by (1301) , so comparison 
of values for different e is problematic. One might argue that the reduction of the sampling area and the related 
decrease in value would warrant a scaling that should depend in one form or the other on e. Ultimately, this kind of 
manipulation is futile, because the A 2 do not depend only on the pupil area but also on its shape.] A noticeable effect 

is the break-up of the degeneracy of the values for and as e increases. It is plausible that modes that do not 
have a node in the middle of the pupil are more sensitive to blanking of the center part than modes that approach 
zero as x — > 0. 




III. SUMMARY 



The text demonstrates computation of the statistically independent modes of phase screens in the von Karman 
class of outer scales sampled over spherical receiver pupils, optionally taking into account masking of a circular region 
in the center. The main value in this work resides in techniques to compute integrals in the eigenvalue problem. 

These modes can be employed to compute outer scales from measured phase functions. The goodness-of-fit of the 
von Karman modal spectrum can be addressed. The number of degrees of freedom depends on how many of the 
modes remain free from other genuine sources of noise besides atmospheric turbulence. 
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FIG. 6: Modal eigenvalues for e — 0.0462 (left, roughly representative of the Keck center tile size) and e = 0.1369 (right) with 
different symbols for inverse outer scales ctq from to 0.8. 



Appendix A: Karhunen-Loeve Modes of the Phase Covariance Inside the Circular Pupil 

Explicit expansions of the first (large-scale) modes (|2"21 follow for dimensionless inverse outer scales from 0.05 to 1. 
The case of infinite outer scale, ctq = 0, has been published before [3l[ and is not duplicated here. Each line shows the 
eigenvalue (|30[) . a colon, the azimuthal parameter another colon, and a sum over Zernike functions using Noll's 

indexing [32]. The sm(q9) and cos(q9) are basically degenerate if q ^ 0: two modes share the same radial (x) 
and Ap and differ only in M q (9). A notation with two comma-separated indices for Z condenses them into a single 

line, and choosing either the left or the right number in the index pairs selects a unique Kp q \x.). 

The squared coefficients sum up to 1/-7T to cancel Noll's normalization of the Zernike basis and to end up with a 
Kp q ^ normalized to unity. 

cr = 0.05: 

0.1270639 : 1: +0 . 5621830*Z(2 , 3) -0 . 0474462*Z(8,7) +0 . 0030025*Z(16 , 17) -0.0000182*Z(30,29) 
+0.0000032*Z(46,47) 

0.0176247 : 2: +0 . 5540843*Z(6 , 5) -0 . 1056709*Z(12 , 13) +0.0115730*Z(24,23) -0.0004187*Z(38,39) 
+0 . 0000243*Z (58 , 57) +0 . 0000014*Z (80 , 81 ) 

0.0176247 : 0: +0 . 5540843*Z(1) -0 . 1056709*Z(4) +0.0115730*Z(11) -0 . 0004187*Z(22) +0.0000243*Z(37) 
+0.0000014*Z(56) 

0.0051337 : 3: +0 . 5411529*Z(10 ,9) -0 . 1576676*Z(18, 19) +0.0245247*Z(32,31) -0 . 0016922*Z(48,49) 
+0.0001077*Z(70,69) 

0.0043151 : 1: +0 . 0457305*Z(2 ,3) +0 . 5301453*Z(8,7) -0 . 1848336*Z(16 , 17) +0.0315406*Z(30,29) 
-0.0025078*Z(46,47) +0 . 0001649*Z(68 ,67) -0 . 0000020*Z(92,93) +0 . 0000011*Z(122 , 121) 

0.0021374 : 4: +0 . 5271594*Z(14, 15) -0. 1971752*Z(26,25) +0.0389824*Z(40,41) -0 . 0038784*Z(60,59) 
+0.0003054*Z(82,83) -0 . 0000094*Z(110, 109) +0 . 0000018*Z(140 , 141) 

0.0016213 : 2: +0 . 0990281*Z(6 , 5) +0 . 4914390*Z(12 , 13) -0.2517952*Z(24,23) +0.0594666*Z(38,39) 
-0.0073175*Z(58,57) +0 . 0006477*Z(80 ,81) -0 . 0000312*Z(108, 107) +0 . 0000032*Z(138 , 139) 



0.0016213 : 0: +0 . 0990281*Z(1) +0.4914390*Z(4) -0.2517952*Z(11) +0 . 0594666*Z(22) -0.0073175*Z(37) 
+0.0006477*Z(56) -0 . 0000312*Z(79) +0.0000032*Z(106) 
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0.0010797 : 5: +0 . 5135000*Z(20 ,21) -0.2273859*Z(34,33) +0.0536290*Z(50,51) -0.0068392*Z(72,71) 
+0.0006548*Z(96,97) -0 . 0000348*Z(126, 125) +0 . 0000037*Z(158 , 159) 

0.0007713 : 3: +0 . 1435201*Z(10 ,9) +0.4445464*Z(18, 19) -0.3025397*Z(32,31) +0.0912672*Z(48,49) 
-0.0150800*Z(70,69) +0.0017271*Z(94,95) -0.0001286*Z(124, 123) +0 . 0000102*Z(156 , 157) 



a = 0.1: 

0.0802920 : 1: +0 . 5603194*Z(2 ,3) -0 . 0658073*Z(8,7) +0 . 0046317*Z(16 , 17) -0.0000506*Z(30,29) 
+0.0000052*Z(46,47) 

0.0155178 : 2: +0 . 5522704*Z(6 , 5) -0 . 1146220*Z(12 , 13) +0.0129920*Z(24,23) -0.0005054*Z(38,39) 
+0 . 0000290*Z (58 , 57) +0 . 0000014*Z (80 , 81 ) 

0.0155178 : 0: +0 . 5522704*Z(1) -0 . 1146220*Z(4) +0.0129920*Z(11) -0 . 0005054*Z(22) +0.0000290*Z(37) 
+0.0000014*Z(56) 

0.0049007 : 3: +0 . 5400346*Z(10 ,9) -0. 1613201*Z(18, 19) +0 . 0254004*Z(32 ,31) -0 . 0017817*Z(48 ,49) 
+0.0001143*Z(70,69) 

0.0039238 : 1: +0 . 0631227*Z(2 ,3) +0 . 5235370*Z(8,7) -0 . 1974575*Z(16 , 17) +0.0351595*Z(30,29) 
-0.0029586*Z(46,47) +0 . 0001991*Z(68 ,67) -0 . 0000034*Z(92,93) +0 . 0000013*Z(122 , 121) 

0.0020876 : 4: +0 . 5264522*Z(14, 15) -0. 1989426*Z(26,25) +0.0395483*Z(40,41) -0 . 0039580*Z(60 ,59) 
+0.0003132*Z(82,83) -0 . 0000098*Z(110, 109) +0 . 0000019*Z(140 , 141) 

0.0015585 : 2: +0 . 1070746*Z(6 , 5) +0 . 4864405*Z(12 , 13) -0.2575283*Z(24,23) +0.0619577*Z(38,39) 
-0.0077766*Z(58,57) +0 . 0006978*Z(80 ,81) -0.0000347*Z(108, 107) +0 . 0000034*Z(138 , 139) 

0.0015585 : 0: +0. 1070746*Z(1) +0 . 4864405*Z(4) -0.2575283*Z(11) +0 . 0619577*Z(22) -0.0077766*Z(37) 
+0.0006978*Z(56) -0 . 0000347*Z(79) +0.0000034*Z(106) 

0.0010641 : 5: +0 . 5130108*Z(20 ,21) -0 . 2283914*Z(34,33) +0.0540274*Z(50,51) -0 . 0069096*Z(72 ,71) 
+0.0006633*Z(96,97) -0 . 0000353*Z(126 , 125) +0 . 0000038*Z(158, 159) 

0.0007576 : 3: +0 . 1466064*Z(10 ,9) +0.4417699*Z(18, 19) -0.3046954*Z(32,31) +0 . 0926064*Z(48,49) 
-0.0154085*Z(70,69) +0 . 0017748*Z(94,95) -0 . 0001332*Z(124, 123) +0 . 0000106*Z(156 , 157) 



(J = 0.15: 

0.0544954 : 1: +0 . 5580974*Z(2 ,3) -0 . 0824405*Z(8,7) +0 . 0063821*Z(16 , 17) -0.0000923*Z(30,29) 
+0.0000081*Z(46,47) 

0.0132897 : 2: +0 . 5500307*Z(6 , 5) -0 . 1247356*Z(12 , 13) +0.0147219*Z(24,23) -0.0006141*Z(38,39) 
+0 . 0000353*Z (58 , 57) +0 . 0000015*Z (80 , 81 ) 

0.0132897 : 0: +0 . 5500307*Z(1) -0 . 1247356*Z(4) +0.0147219*Z(11) -0 . 0006141*Z(22) +0.0000353*Z(37) 
+0.0000015*Z(56) 

0.0045724 : 3: +0 . 5384109*Z(10 ,9) -0. 1664591*Z(18, 19) +0 . 0266702*Z(32 ,31) -0 . 0019125*Z(48,49) 
+0 . 0001243*Z (70 , 69) +0 . 0000010*Z(124 , 123) 

0.0035533 : 1: +0 . 0787021*Z(2 ,3) +0 . 5165478*Z(8,7) -0 . 2092330*Z(16 , 17) +0.0387812*Z(30,29) 
-0.0034279*Z(46,47) +0 . 0002365*Z(68,67) -0 . 0000049*Z(92,93) +0 . 0000014*Z(122 , 121) 

0.0020105 : 4: +0 . 5253405*Z(14, 15) -0 . 2016812*Z(26,25) +0.0404365*Z(40,41) -0 . 0040835*Z(60 , 59) 
+0.0003256*Z(82,83) -0 . 0000104*Z(110, 109) +0 . 0000019*Z(140 , 141) 

0.0014813 : 2: +0. 1160947*Z(6,5) +0 . 4804411*Z(12 , 13) -0.2640555*Z(24,23) +0.0649074*Z(38,39) 
-0.0083296*Z(58,57) +0 . 0007596*Z(80 ,81) -0 . 0000389*Z(108, 107) +0 . 0000037*Z(138 , 139) 

0.0014813 : 0: +0 . 1160947*Z(1) +0 . 4804411*Z(4) -0.2640555*Z(11) +0 . 0649074*Z(22) -0.0083296*Z(37) 
+0.0007596*Z(56) -0 . 0000389*Z(79) +0.0000037*Z(106) 
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0.0010391 : 5: +0 . 5122189*Z(20 ,21) -0.2300063*Z(34,33) +0.0546713*Z(50,51) -0.0070238*Z(72,71) 
+0.0006773*Z(96,97) -0.0000363*Z(126, 125) +0 . 0000039*Z(158 , 159) 

0.0007373 : 3: +0. 1509283*Z(10,9) +0.4377755*Z(18, 19) -0.3077212*Z(32,31) +0.0945258*Z(48,49) 
-0.0158828*Z(70,69) +0.0018443*Z(94,95) -0.0001398*Z(124, 123) +0 . 0000111*Z(156 , 157) 



cr = 0.2: 

0.0384868 : 1: +0 . 5555830*Z(2 ,3) -0 . 0978254*Z(8,7) +0 . 0082237*Z(16 , 17) -0.0001412*Z(30,29) 
+0 . 00001 18*Z (46 ,47) +0 . 000001 1*Z (68 , 67) 

0.0112251 : 2: +0 . 5475076*Z(6 , 5) -0 . 1351598*Z(12 , 13) +0.0166326*Z(24,23) -0.0007383*Z(38,39) 
+0 . 0000432*Z (58 , 57) +0 . 0000017*Z (80 , 81 ) 

0.0112251 : 0: +0 . 5475076*Z(1) -0 . 1351598*Z(4) +0.0166326*Z(11) -0 . 0007383*Z(22) +0.0000432*Z(37) 
+0.0000017*Z(56) 

0.0041940 : 3: +0 . 5364214*Z(10 ,9) -0. 1725149*Z(18, 19) +0.0282171*Z(32,31) -0 . 0020739*Z(48,49) 
+0.0001371*Z(70,69) +0 . 0000011*Z(124, 123) 

0.0032042 : 1: +0 . 0929350*Z(2 ,3) +0 . 5092242*Z(8,7) -0 . 2203375*Z(16 , 17) +0.0424228*Z(30,29) 
-0.0039175*Z(46,47) +0 . 0002775*Z(68,67) -0 . 0000066*Z(92,93) +0 . 0000016*Z(122 , 121) 

0.0019126 : 4: +0 . 5238905*Z(14, 15) -0 . 2051844*Z(26,25) +0.0415913*Z(40,41) -0 . 0042476*Z(60 ,59) 
+0.0003420*Z(82,83) -0 . 0000111*Z(110 , 109) +0 . 0000020*Z(140 , 141) 

0.0013966 : 2: +0. 1252979*Z(6,5) +0 . 4738551*Z(12 , 13) -0.2708383*Z(24,23) +0.0681004*Z(38,39) 
-0.0089399*Z(58,57) +0 . 0008295*Z(80 ,81) -0 . 0000438*Z(108, 107) +0 . 0000041*Z(138 , 139) 

0.0013966 : 0: +0. 1252979*Z(1) +0 . 4738551*Z(4) -0.2708383*Z(11) +0 . 0681004*Z(22) -0.0089399*Z(37) 
+0.0008295*Z(56) -0 . 0000438*Z(79) +0.0000041*Z(106) 

0.0010062 : 5: +0 . 5111531*Z(20 ,21) -0 . 2321555*Z(34,33) +0.0555358*Z(50,51) -0 . 0071778*Z(72 ,71) 
+0.0006962*Z(96,97) -0 . 0000375*Z(126 , 125) +0 . 0000040*Z(158, 159) 

0.0007122 : 3: +0 . 1559847*Z(10 ,9) +0.4329305*Z(18, 19) -0.3112822*Z(32,31) +0 . 0968427*Z(48,49) 
-0.0164612*Z(70,69) +0 . 0019301*Z(94,95) -0 . 0001480*Z(124, 123) +0 . 0000118*Z(156 , 157) 



cr = 0.3: 

0.0207627 : 1: +0 . 5499437*Z(2 ,3) -0. 1254079*Z(8,7) +0 . 0120252*Z(16 , 17) -0.0002530*Z(30,29) 
+0 . 0000226*Z (46 ,47) +0 . 0000019*Z (68 , 67) 

0.0078797 : 2: +0 . 5419000*Z(6 , 5) -0 . 1556381*Z(12 , 13) +0.0207357*Z(24,23) -0.0010191*Z(38,39) 
+0 . 0000634*Z (58 , 57) +0 . 0000021*Z (80 , 81) 

0.0078797 : 0: +0 . 5419000*Z(1) -0 . 1556381*Z(4) +0.0207357*Z(11) -0 . 0010191*Z(22) +0.0000634*Z(37) 
+0.0000021*Z(56) 

0.0034105 : 3: +0 . 5317006*Z(10 ,9) -0. 1859701*Z(18, 19) +0 . 0318335*Z(32 ,31) -0 . 0024607*Z(48,49) 
+0.0001696*Z(70,69) -0 . 0000011*Z(94,95) +0 . 0000013*Z(124, 123) 

0.0025805 : 1: +0 . 1179237*Z(2 ,3) +0 . 4939434*Z(8,7) -0 . 2406863*Z(16 , 17) +0.0496921*Z(30,29) 
-0.0049451*Z(46,47) +0 . 0003693*Z(68,67) -0 . 0000106*Z(92,93) +0 . 0000020*Z(122 , 121) 

0.0016809 : 4: +0 . 5202115*Z(14, 15) -0 . 2137474*Z(26,25) +0.0444970*Z(40,41) -0.0046666*Z(60,59) 
+0.0003854*Z(82,83) -0 . 0000130*Z(110 , 109) +0 . 0000023*Z(140 , 141) 

0.0012199 : 2: +0 . 1430463*Z(6 , 5) +0 . 4596855*Z(12 , 13) -0.2843055*Z(24,23) +0.0748388*Z(38,39) 
-0.0102679*Z(58,57) +0 . 0009874*Z(80 ,81) -0.0000551*Z(108, 107) +0 . 0000049*Z(138 , 139) 

0.0012199 : 0: +0 . 1430463*Z(1) +0 . 4596855*Z(4) -0.2843055*Z(11) +0 . 0748388*Z(22) -0.0102679*Z(37) 
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+0.0009874*Z(56) -0 . 0000551*Z(79) +0.0000049*Z(106) 

0.0009225 : 5: +0 . 5083277*Z(20 ,21) -0.2377231*Z(34,33) +0.0578144*Z(50,51) -0.0075872*Z(72,71) 
+0.0007473*Z(96,97) -0.0000409*Z(126, 125) +0 . 0000043*Z(158 , 159) 

0.0006527 : 3: +0. 1670529*Z(10,9) +0.4216099*Z(18, 19) -0.3191848*Z(32,31) +0. 1022108*Z(48,49) 
-0.0178264*Z(70,69) +0.0021369*Z(94,95) -0.0001680*Z(124, 123) +0 . 0000134*Z(156 , 157) 



a = 0.4: 

0.0120802 : 1: +0 . 5438669*Z(2 ,3) -0 . 1492296*Z(8,7) +0 . 0157816*Z(16 , 17) -0.0003729*Z(30,29) 
+0 . 0000384*Z (46 , 47) +0 . 0000031*Z (68 , 67) 

0.0055228 : 2: +0 . 5358544*Z(6 , 5) -0 . 1747716*Z(12 , 13) +0.0249608*Z(24,23) -0.0013265*Z(38,39) 
+0.0000892*Z(58,57) +0 . 0000027*Z(80,81) +0 . 0000011*Z(108, 107) 

0.0055228 : 0: +0 . 5358544*Z(1) -0 . 1747716*Z(4) +0.0249608*Z(11) -0 . 0013265*Z(22) +0.0000892*Z(37) 
+0 . 0000027*Z (56) +0 . 0000011*Z (79) 

0.0027015 : 3: +0 . 5263645*Z(10 ,9) -0. 1998967*Z(18, 19) +0 . 0358207*Z(32 ,31) -0 . 0029028*Z(48 ,49) 
+0.0002098*Z(70,69) -0 . 0000017*Z(94,95) +0 . 0000016*Z(124, 123) 

0.0020613 : 1: +0 . 1388056*Z(2 ,3) +0.4783720*Z(8,7) -0 . 2587251*Z(16 , 17) +0.0568160*Z(30,29) 
-0.0060141*Z(46,47) +0.0004724*Z(68,67) -0.0000151*Z(92,93) +0 . 0000025*Z(122 , 121) 

0.0014359 : 4: +0 . 5158045*Z(14, 15) -0 . 2234543*Z(26,25) +0.0479255*Z(40,41) -0 . 0051719*Z(60 ,59) 
+0.0004400*Z(82,83) -0 . 0000153*Z(110 , 109) +0 . 0000027*Z(140 , 141) 

0.0010481 : 2: +0. 1591466*Z(6,5) +0 . 4448988*Z(12 , 13) -0.2970362*Z(24,23) +0.0817255*Z(38,39) 
-0.0116808*Z(58,57) +0 . 0011637*Z(80 ,81) -0 . 0000680*Z(108, 107) +0 . 0000059*Z(138 , 139) 

0.0010481 : 0: +0. 1591466*Z(1) +0 . 4448988*Z(4) -0.2970362*Z(11) +0 . 0817255*Z(22) -0.0116808*Z(37) 
+0.0011637*Z(56) -0 . 0000680*Z(79) +0.0000059*Z(106) 

0.0008254 : 5: +0 . 5047821*Z(20 ,21) -0 . 2444603*Z(34,33) +0.0606444*Z(50,51) -0.0081031*Z(72,71) 
+0.0008134*Z(96,97) -0 . 0000452*Z(126, 125) +0 . 0000048*Z(158, 159) 

0.0005874 : 3: +0 . 1782251*Z(10 ,9) +0.4090493*Z(18, 19) -0.3273494*Z(32,31) +0 . 1081019*Z(48 ,49) 
-0.0193653*Z(70,69) +0 . 0023770*Z(94,95) -0 . 0001917*Z(124, 123) +0 . 0000154*Z(156 , 157) 



(T = 0.6: 

0.0048075 : 1: +0 . 5316377*Z(2 ,3) -0 . 1875067*Z(8,7) +0 . 0226296*Z(16 , 17) -0.0006127*Z(30,29) 
+0.0000863*Z(46,47) +0.0000078*Z(68,67) +0 . 0000016*Z(92 ,93) 

0.0028053 : 2: +0 . 5234235*Z(6 , 5) -0 . 2079357*Z(12 , 13) +0.0331139*Z(24,23) -0.0019680*Z(38,39) 
+0.0001565*Z(58,57) +0 . 0000053*Z(80,81) +0 . 0000020*Z(108, 107) 

0.0028053 : 0: +0 . 5234235*Z(1) -0 . 2079357*Z(4) +0.0331139*Z(11) -0 . 0019680*Z(22) +0.0001565*Z(37) 
+0 . 0000053*Z (56) +0 . 0000020*Z (79) 

0.0016497 : 3: +0 . 5149116*Z(10 ,9) -0.2263237*Z(18, 19) +0 . 0440286*Z(32 ,31) -0 . 0038640*Z(48 ,49) 
+0.0003090*Z(70,69) -0 . 0000025*Z(94,95) +0 . 0000025*Z(124, 123) 

0.0013057 : 1: +0 . 1704452*Z(2 ,3) +0 . 4481773*Z(8,7) -0 . 2887852*Z(16 , 17) +0.0702216*Z(30,29) 
-0.0081810*Z(46,47) +0 . 0007039*Z(68 ,67) -0 . 0000250*Z(92,93) +0 . 0000039*Z(122 , 121) 

0.0009994 : 4: +0 . 5058217*Z(14, 15) -0 . 2435913*Z(26,25) +0.0554728*Z(40,41) -0 . 0063270*Z(60 ,59) 
+0.0005743*Z(82,83) -0 . 0000206*Z(110 , 109) +0 . 0000037*Z(140 , 141) 

0.0007519 : 2: +0 . 1855962*Z(6 , 5) +0 . 4154116*Z(12 , 13) -0.3194096*Z(24,23) +0.0951789*Z(38,39) 
-0.0146004*Z(58,57) +0 . 0015527*Z(80 ,81) -0 . 0000970*Z(108, 107) +0 . 0000085*Z(138 , 139) 
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0.0007519 : 0: +0. 1855962*Z(1) +0.4154116*Z(4) -0.3194096*Z(11) +0 . 0951789*Z(22) -0 . 0146004*Z(37) 
+0.0015527*Z(56) -0 . 0000970*Z(79) +0.0000085*Z(106) 

0.0006290 : 5: +0 . 4963318*Z(20 ,21) -0.2595260*Z(34,33) +0.0672512*Z(50,51) -0.0093399*Z(72,71) 
+0.0009791*Z(96,97) -0.0000558*Z(126, 125) +0 . 0000061*Z(158 , 159) 

0.0004592 : 3: +0 . 1984313*Z(10 ,9) +0.3826865*Z(18, 19) -0.3427694*Z(32,31) +0. 1202986*Z(48,49) 
-0.0226871*Z(70,69) +0.0029188*Z(94,95) -0.0002463*Z(124, 123) +0 . 0000203*Z(156 , 157) 



a = 0.8: 

0.0022407 : 1: +0 . 5203650*Z(2 ,3) -0 . 2161624*Z(8,7) +0 . 0283406*Z(16 , 17) -0.0008481*Z(30,29) 
+0.0001528*Z(46,47) +0 . 0000171*Z(68 ,67) +0 . 0000034*Z(92 ,93) 

0.0015169 : 2: +0 . 5114918*Z(6 , 5) -0 . 2346225*Z(12 , 13) +0.0403918*Z(24,23) -0.0025985*Z(38,39) 
+0.0002410*Z(58,57) +0 . 0000108*Z(80 ,81) +0 . 0000035*Z(108, 107) 

0.0015169 : 0: +0 . 5114918*Z(1) -0 . 2346225*Z(4) +0.0403918*Z(11) -0 . 0025985*Z(22) +0.0002410*Z(37) 
+0 . 0000108*Z (56) +0 . 0000035*Z (79) 

0.0010125 : 3: +0 . 5034249*Z(10 ,9) -0.2493257*Z(18, 19) +0 . 0518292*Z(32 ,31) -0 . 0048426*Z(48,49) 
+0.0004262*Z(70,69) -0 . 0000017*Z(94,95) +0 . 0000039*Z(124, 123) 

0.0008339 : 1: +0 . 1919391*Z(2 ,3) +0 . 4207723*Z(8,7) -0 . 3123424*Z(16 , 17) +0.0821921*Z(30,29) 
-0.0102809*Z(46,47) +0 . 0009559*Z(68 ,67) -0 . 0000348*Z(92,93) +0 . 0000058*Z(122 , 121) 
+0.0000011*Z(154,155) 

0.0006790 : 4: +0 . 4953154*Z(14, 15) -0 . 2625515*Z(26,25) +0.0631000*Z(40,41) -0 . 0075552*Z(60 ,59) 
+0.0007312*Z(82,83) -0 . 0000257*Z(110 , 109) +0 . 0000052*Z(140 , 141) +0 . 0000010*Z(176, 175) 

0.0005307 : 2: +0 . 2050703*Z(6 , 5) +0 . 3877207*Z(12 , 13) -0.3376902*Z(24,23) +0. 1076308*Z(38,39) 
-0.0174870*Z(58,57) +0 . 0019688*Z(80 ,81) -0 . 0001283*Z(108, 107) +0 . 0000118*Z(138, 139) 

0.0005307 : 0: +0 . 2050703*Z(1) +0 . 3877207*Z(4) -0.3376902*Z(11) +0 . 1076308*Z(22) -0.0174870*Z(37) 
+0.0019688*Z(56) -0 . 0001283*Z(79) +0.0000118*Z(106) 

0.0004624 : 5: +0 . 4870068*Z(20,21) -0 . 2747635*Z(34,33) +0.0743178*Z(50,51) -0 . 0107141*Z(72 ,71) 
+0.0011749*Z(96,97) -0 . 0000677*Z(126, 125) +0 . 0000079*Z(158, 159) 

0.0003491 : 3: +0 . 2146554*Z(10 ,9) +0.3567923*Z(18, 19) -0.3560623*Z(32,31) +0 . 1321324*Z(48,49) 
-0.0260852*Z(70,69) +0 . 0035049*Z(94,95) -0 . 0003067*Z(124, 123) +0 . 0000262*Z(156 , 157) 



a = 1.0: 

0.0011717 : 1: +0 . 5104798*Z(2 ,3) -0 . 2379706*Z(8,7) +0 . 0329996*Z(16 , 17) -0.0010862*Z(30,29) 
+0.0002291*Z(46,47) +0 . 0000318*Z(68 ,67) +0 . 0000066*Z(92,93) +0 . 0000016*Z(122 , 121) 

0.0008736 : 2: +0 . 5006014*Z(6 , 5) -0 . 2559705*Z(12 , 13) +0.0466574*Z(24,23) -0.0031965*Z(38,39) 
+0.0003360*Z(58,57) +0 . 0000203*Z(80 ,81) +0 . 0000060*Z(108, 107) +0 . 0000015*Z(138, 139) 

0.0008736 : 0: +0 . 5006014*Z(1) -0 . 2559705*Z(4) +0.0466574*Z(11) -0 . 0031965*Z(22) +0.0003360*Z(37) 
+0.0000203*Z(56) +0 . 0000060*Z(79) +0.0000015*Z(106) 

0.0006374 : 3: +0 . 4925470*Z(10 ,9) -0.2687151*Z(18, 19) +0 . 0588705*Z(32 ,31) -0 . 0057865*Z(48 ,49) 
+0.0005542*Z(70,69) +0 . 0000017*Z(94,95) +0.0000061*Z(124, 123) +0 . 0000015*Z(156 , 157) 

0.0005433 : 1: +0 . 2065358*Z(2 ,3) +0 . 3967198*Z(8,7) -0 . 3309519*Z(16 , 17) +0.0926588*Z(30,29) 
-0.0122456*Z(46,47) +0.0012152*Z(68,67) -0 . 0000426*Z(92,93) +0 . 0000086*Z(122 , 121) 
+0.0000016*Z(154,155) 

0.0004619 : 4: +0 . 4849996*Z(14, 15) -0 . 2794018*Z(26 ,25) +0.0703001*Z(40,41) -0 . 0087740*Z(60 , 59) 
+0.0009009*Z(82,83) -0 . 0000295*Z(110, 109) +0 . 0000073*Z(140 , 141) +0 . 0000014*Z(176, 175) 
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0.0003744 : 2: +0 . 2190750*Z(6 , 5) +0 . 3626864*Z(12 , 13) -0.3524757*Z(24,23) +0. 1188288*Z(38,39) 
-0.0202349*Z(58,57) +0.0023925*Z(80,81) -0 . 0001595*Z(108, 107) +0 . 0000158*Z(138 , 139) 
+0.0000012*Z(174,173) 

0.0003744 : 0: +0 . 2190750*Z(1) +0.3626864*Z(4) -0.3524757*Z(11) +0. 1188288*Z(22) -0 . 0202349*Z(37) 
+0.0023925*Z(56) -0 . 0001595*Z(79) +0.0000158*Z(106) +0 . 0000012*Z(137) 

0.0003356 : 5: +0.4775241*Z(20,21) -0.2890089*Z(34,33) +0.0812749*Z(50,51) -0.0121207*Z(72,71) 
+0.0013880*Z(96,97) -0 . 0000795*Z(126, 125) +0 . 0000103*Z(158 , 159) +0 . 0000013*Z(196, 195) 

0.0002623 : 3: -0 . 2270448*Z(10 ,9) -0.3325947*Z(18, 19) +0.3671201*Z(32,31) -0. 1431189*Z(48,49) 
+0.0293962*Z(70,69) -0.0041058*Z(94,95) +0.0003695*Z(124, 123) -0 . 0000330*Z(156 , 157) 



Appendix B: Piston Removal 



The piston mode over the annulus is a mode without azimuth dependence, 5 = 0, where K p °\x) is a constant that 
does not depend on the radial coordinate x. Finite outer scales, k Q > 0, prohibit the ultraviolet divergence of the 
integrals that is characteristic to the Kolmogorov variance, so piston subtraction could be implemented by ignoring 
the associated mode. Yet, removal of the piston improves numerical stability because this contribution otherwise 
dominates the matrix elements at q = 0. 

At q = 0, equation (f3"8j) reads 



k 



Jo(kx)J (kx') = ^(-l)"(2n + l)J 2n+1 (fe) 



n>0 
X 2*1 



-n, n + 1 
1 



2*1 



-n, n + 1 , . 9 ^ 
' 1 | sin $ 



(Bl) 



The dependence on x and x' is generated by the product of the 2*1 (•), which is a symmetric bi-variate polynomial in 
x and x' of mixed order 2n. The basic examples are 



2*1 



-n, n + 1 
1 



2(x 2 

f„2 , 



2*1 



- X 
„'2 



-n, n + 1 
1 



4x 2 x' 2 ) 



sin 2 $ 



n=0, 
n=l, 
n=2. 



(B2) 



-1 + 12(x 2 + x 12 ) - 30(x 4 + x' 4 + 4x 2 x' 2 ) + 20(x 6 + x m + §x A x' 2 + §x 2 x' 4 ) n=3. 



The simplicity of this format might be unexpected in view of the square roots in (|4ip . but if ()38j) is constructed 
term-by-term on the right hand side by balancing the powers in k on both sides of the Neuman series, the linear 
coupling of k and x, x' in the arguments on the left hand side makes this plausible. 

The piston contribution is manipulated by adding some 8 n to these products, which by Mercer's theorem effects 
the piston's X 2 . To suppress the piston, equation (|57)l and the orthogonality of the Kp°\x) suggest that 



'dx' 



udx 



2^1 



-n, n + 1 
1 



sin 



2-fl 



-n, n + 1 
1 



sin 2 $ + 6 n 



= 0. 



(B3) 



Insertion of {B2| yields for example 



r -1 

-l-2e 2 

-5(l + 2e 2 + 2e 4 ) 
-31(1 + 2e 2 ) - 10e 4 ' 



+ 7e 2 ) 



-229(1 + 2e 2 ) - 2e 4 (334 + 329e 2 + 294e 4 ) 
I -1891(1 + 2e 2 ) - 14e 4 (403 + 488e 2 + 408e 4 



396e 6 



n=0, 
n=l, 
n=2, 
n=3, 
n=4, 
n=5. 



(B4) 
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The first value removes the entire term at n — in (|B1|) . The results in the main section added the 8 n to all n < 10. 
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